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We estimate the energy cost associated with two pancake vortices colliding in a layered supercon- 
l/") . ductor. It is argued that this energy sets the plastics energy scale and is the analogue of the crossing 

0^ ' energy for vortices in the continuum case. The starting point of the calculation is the Lawrence- 

Doniach version of the Ginzburg-Landau free energy for type-II superconductors. The magnetic 
fields considered are along the c-direction and assumed to be sufficiently high that the lowest Lan- 
dau level approximation is valid. For Bi-2212, where it is know that layering is very important, 
the results are radically different from what would have been obtained using a three-dimensional 
anisotropic continuum model. We then use the plastic energy for Bi-2212 to successfully explain 
recent results from Hellerqvist et al. on its longitudinal resistance. 
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There has beeOpHiuch recent experimental interest in the dependence of the longitudinal resistivity in the HTC- 
superconductorsata on field and temperature. In a recent paper, we produced a theory for some aspects of this data 
which involved knowing the crossing energy for two vortices (flux lines) . Our work was essentially an extension of the 
(T) . hydrodynamic_theory of Marchetti, Nelson and CatesSi3. It gave a reasonable fit to the longitudinal resistivity data 
of Safar et alu on YBCO crystals. However, it was clear from our calculation that the length scale over which the 
crossing took place was of the order of the spacing between the superconducting Cu-0 planes. Our calculation had 
been done using anisotropic continuum Ginzburg-Landau theory, and so neglected any layering effects. The crossing 
energy that we calculated for the three-dimensional anisotropic case can be thought of as the plastic energy cost of 

S making non-trivial topological deformations of the vortices, and we shall refer to it as U^f nt . 
In this paper we calculate the energy associated with two pancakes colliding within a layer: this is an estimate of the 
. plastic deformation energy for vortices in a strongly layered superconductor and sets the timescales for deformation 
of the vortices within the system. Intuitively one might have expected that neighbouring pancake vortices in a layer 
would have exchanged positions by moving around each other. However, within the Lowest Landau Level (LLL) 
approximation which we use, Dodgson and MooreQ have shown that the lowest energy barrier for this to occur is when 
, the vortices collide and pass through each other. 

We start from a Lawrence-Doniach Ginzburg-Landau free energy functional We find that including the layering 
leaves the results for the plastic energy of YBCO essentially unchanged but gives markedly different results for Bi-2212. 
This suggests that treating YBCO as a 3D continuum anisotropic superconductor within GL theory is reasonable, 
but that in Bi-2212 the effect of the Cu-0 layers must be included — as was to be expected as the coupling between 
the layers in Bi-2212 is known to be very weak. 

To fix notations, we briefly describe the modified Ginzburg-Landau theory for a layered superconductor, where ip, 
is our spatially dependent order parameter. We begin with the continuum free energy functional, 




Here a(T) is the temperature-dependent variable, f3 is the coupling constant, and are the effective masses. In 
the cases we consider the masses in the afr-plane are taken as equal and are denoted by m a b, and the mass in the 
c-direction is written as m c . The temperature dependence of a(T) is taken to be of the form , a(T) = (T — T c )a' . 
We also assume the LLL approximation (see Refs. |9|Jl~0| for justification of it). In this approximation the derivative 
terms in the direction perpendicular to the planes reduce to ehB /m a b\4>\ 2 , so it is convenient to write the temperature 
variable as an a + eBh/m a b- This is zero along the H C2 line, negative below and positive above. According to 
Tesanovic et al.EiJ the LLL approximation can be trusted when H > H C2 / f, where / = 3. 

In our previous calculation we chose to consider a system that contained only two vortices, but with different types 
of periodic boundary conditions. In this paper we instead generalise to the layered case the procedure of Dodgson 
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and Moorda who allowed only the two vortices that are crossing to move and kept all the other vortices fixed in 
their triangular lattice positions. The advantage of using a background of fixed vortices is that the resulting crossing 
energy, U p \ is an upper-bound on the true crossing energy. Furthermore, it has been shown that the corrections to 
this bound seem to be very small, at least for the continuum limitB. 

The general form for the order parameter in the afr-plane within the LLL approximation is: 

4>(x,y) = y (C)e _ ^ IC|2jB/C2 * o) where £ = x + iy, (2) 
and /(C) can be written in the product form, 

/(OcxII^C-O) (3) 

where the Ci are the (complex) positions of the vortices. The minimum potential free energy configuration for this 
system is the triangular lattice of lattice spacing -\/3/2Z 2 = <&o/B. We will label the order-parameter that describes 
this configuration by ipo(x,y). 

We are interested in two vortices being displaced along the x-axis from their lattice positions at £ = ±1/2 and 
colliding at the origin x = y = 0. Using the formalism of Ref. || we have for the order parameter at height z 

^ w) = M X J C+ pl%* {z)) (4) 

2a(z) is the separation of the colliding vortices and a(0) = 0. (see Fig. |l|). 

Layering will be taken into account by replacing in Eq. [l] the derivative dip(x, y, z)/dz by (ip n +i(x, y) —ip n (x,y))/d 
where ip n (x,y) denotes the value of ip in the n th layer at point (x,y) and d is the effective spacing between the 
superconducting Cu-0 layers. In the context of the LLL formalism, such an approximation has been derived as a 
consequence of a periodically varying T c (z) by Sasik and Stroudli3. Note that within the LLL approximation the 
vector potential A is always that of the uniform field so A z = 

The resulting effective Ginzburg-Landau free energy functional within the LLL approximation for a layered structure 
can then be written as, 

F = d T,j (-l^ll^l' + flV^J' + ^^l^+i-^l^dxdy (5) 

The previous results for the continuum case of Dodgson and Moorei, with a rigid lattice background of vortices led to 
a crossing energy of 2.3 h&o \ceii\ 3 ^ 2 /{\J2m c f3B), which is a rigorous upper bound to the crossing energy. Our earlier cal- 
culations using only two vortices with periodic boundary conditions led to an energy of 1.46 ?j$o \otn | 3 ^ 2 / (y/2m c /3B), 
which is somewhat lower, possibly due to errors from the finite size of the system. In this paper we use the rigid 
lattice approach, as the equations are more readily physically interpreted, and easier to solve. 

The extra free energy, AF associated with the two displaced vortices as compared witfuthe original lattice is found 
by substituting Eq. (|4|) into Eq. (||) and then doing the integral over the (x,y) directionsQ, 

= 9(V) where g(fj) = > ] \ > ] c m a^ ] + ^(a 2 n+1 - c? n f (6) 

where a n is the value of a(z) in the n th layer of which there are (2N + 1). The factor (3a — 1.16 is the Abrikosov 
factor. The quantity r\ — d/£ c where £ c = ft/ \/2m c \oLii\ — the 'c-axis correlation length' measures the importance 
of layering effects. The formula for £ c in terms of an should be modifiedtl very close to the H C2 line a regime which 
we therefore shall not consider. When 77 is large, layering effects are important; rj — > is the continuum limit. 
The constants c m and I\ are determined by integration in the xy-plane, cq — 0.7903, C2 = —5.201, C3 = —32.7365, 
C4 = 100.8232, and I\ = 23.5502. The a and (3 terms can be respectively thought of as the 'potential' and 'kinetic' 
contributions to the plastic energy. 

In order for the vortices to collide we require that at the central layer (n = 0) the displaced vortices lie on top of each 
other. By symmetry this should occur at ao = 0. (See Fig. Q.) Our other boundary condition is that a± OQ = ±1/2, 
that is we require the vortices to return to their equilibrium positions in the original lattice. 

The energy barrier for vortices to collide is then found by minimising Eq. (|^), with the boundary conditions imposed. 
The results unlike in the continuum case cannot be written in terms of a reduced temperature variable. They are now 
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field, temperature and layer thickness d dependent. The easiest way to get an idea of how the plastic energy in the 
layered case varies is to compare it with the plastic/crossing energy calculated for the rigid lattice background in the 
continuum case. 

^ydisc 

TT^t = O^fa) (7) 

where C/pi' sc is the crossing energy found in the discrete case and £/p° nt in the continuum case. The form of this 
function is shown in Fig. ||. Putting in typical system parameters for YBCO: \dB/dT\ = 2T/K, m a b/rn c = e 2 = 1/59, 
T c = 93K, and a spacing between Cu-0 planes of 11. 4A (where we have treated the closely spaced pair of Cu-0 layers 
as a 'single' layer) we find that r, ~ 2, for values of field and temperature in the vicinity of the irreversibility line. As 
can be seen from Fig. [| for rj ~ 2 the plastic energy calculated according to the continuum approximation and that 
calculated allowing for layering are approximately the same, indicating that for YBCO the continuum approximation 
is satisfactory, certainly for all the regions of the phase diagram in which we would expect the LLL approximation 
to be valid. The material for which layering effects would be expected to be more important is Bi-2212, where it is 
known from experiments that the coupling between the Cu-0 layers is very weak. Our calculation bears this out; 
using system parameters of |cLB/dT| = 0.5T/K, T c = 86, e = 1/55 (there appears to be a factor of 1/3 uncertainty in 
the value of e!) and d = 30. 9A, we find rj ~ 25 and that the [/distyjycont ^ g_ Thus the continuum estimate appears 
to seriously underestimate the plastic energy for this weakly coupled layered system. 

Analysis of g(rf) shows that for the values of r\ appropriate to Bi-2212 we can greatly simplify the formula for the 
plastic energy, Eq. (^). The displacement of the pancake vortices for these values of 77, is completed within three 
layers (see Fig. ||), and is such that the 'kinetic energy' term in Eq. (||) is negligible in comparison to the 'potential 
energy' term, which tends to the constant cq. Hence when 77 is large 
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In the three dimensional continuum case we were able to use the formula of Cates et al. which says that the 
length scale over which the system is coherent along the field direction, S c can be written as 

6 D ~l e exp[Ug at /KT\. (9) 

The distance l e is the distance along a flux line that is traveled before encountering another flux line, and is related to 
the equilibrium flux line spacing by l e — clq/I p , where l p is a persistence lengthQ which in our case is equal to the c-axis 
correlation length £ c . (There is uncertainty in the numerical constant multiplying l e but it should be of 0(1)). We 



motivated Eq. (|9j) in Ref. 13 in some detail but it can also be readily understood using an explanation of Marchetti 
and NelsonEJ. If we concentrate on the dynamics of one flux line in a flux liquid and label the c-direction by time 
then l/l e becomes an effective attempt frequency to cross or collide, which occurs with probability exp[— U^f nt /kT], 
leading immediately to Eq. (||). Crossing of flux lines always produces a loss of phase coherence. In the layered case 
we can no longer consider 'crossing energies'. However, the arguments detailed above follow through once we consider 
the plastic deformation energy rather than the crossing energy. The importance of the plastic energy for flux line 
motion in the liquid state is shown in Blatter et al&B. n 

We are now in a position to compare with the longitudinal resistivity measurements of Hellerqvist et alu on Bi-2212. 
Using a single crystal of size 2 x 2 x 0.030mm (smallest dimension in the c-direction) they measured the dynamic 
resistance dV/dl as a function of bias current in a magnetic field along the c-direction. Their data of p c against T 
has a smooth peaked behavior with the position of the peak shifting to higher temperatures as the field is lowered. 
By fitting to a power law of the form dV/dl cx I a they find linear I — V behavior for T > T pea k but for T < T pea k 
the exponent a is positive and varies smoothly with temperature and field. 

The apparent disappearance of the longitudinal linear resistivity at the peak suggest that the sample has become 
phase coherent along the c-direction. To locate the point at which the longitudinal resistivity vanishes we set 5 C to 
equal the system size, L = 0.03mm. 

Analysis of the Hellerqvist et ala data, shows that the field-temperature dependence of the longitudinal resistivity 
peaks is well fitted by the formula L ~ S c (see Fig. ||) . At high fields and low temperatures we begin to see deviations 
from the experimental data — but this is no surprise, as these points are so far away from the H C2 line that it is 
unlikely that the LLL approximation remains valid. 

We conclude by observing that the onset of vanishing longitudinal resistance is well-explained by equating the 
plastic energy length scale S c to the system size. No phase transitior£3 need be invoked to explain the experimental 
data. 
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FIG. 1. Trajectory of two lines colliding in the rcy-planc 



FIG. 2. The ratio Up} sc /Up° nt versus r\ = d/£ c - As rj — > the ratio tends to unity while it is linear for large rj. 



FIG. 3. Configuration for vortices colliding in Bi-2212. The positions of the vortices are only defined in the layers 



FIG. 4. Comparison of experiment (points) and theory (solid) line for the line marking the disappearance of the longitudinal 
resistivity in Bi-2212. 
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• Fixed vortices 

O Equilibrium positions of crossing vortices 
O Displaced vortices 
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